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Key Points: 

• Correlated k-distribution models for the VIIRS bands are developed 

• A fast VIIRS instrument simulator is developed 

• A case study is performed to validate the simulator 

Abstract 

A fast instrument simulator is developed to simulate the observations made in cloudy 
atmospheres by the Visible Infrared Imaging Radiometer Suite (VIIRS). The correlated k- 
distribution (CKD) technique is used to compute the transmissivity of absorbing 
atmospheric gases. The bulk scattering properties of ice clouds used in this study are based 
on the ice model used for the MODIS Collection 6 ice cloud products. Two fast radiative 
transfer models based on pre-computed ice cloud look-up-tables are used for the VIIRS 
solar and infrared channels. The accuracy and efficiency of the fast simulator are quantify 
in comparison with a combination of the rigorous line-by-line (LBLRTM) and discrete 
ordinate radiative transfer (DISORT) models. Relative errors are less than 2% for simulated 
TO A reflectances for the solar channels and the brightness temperature differences for the 
infrared channels are less than 0.2 K. The simulator is over three orders of magnitude faster 
than the benchmark LBLRTM+DISORT model. Furthermore, the cloudy atmosphere 
reflectances and brightness temperatures from the fast VIIRS simulator compare favorably 
with those from VIIRS observations. 
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1. Introduction 


The Visible Infrared Imaging Radiometer Suite (VIIRS) on board the Suomi National 
Polar-orbiting Partnership (Suomi NPP) satellite provides critical data for accurately 
determining cloud and aerosol properties, ocean color, sea and land surface temperatures, 
ice motion and temperature, fires, and Earth’s albedo [Lewis et ah, 2010; Lee et ah, 2010; 
Hillger et ah, 2013]. VIIRS includes 16 moderate-resolution channels (referred to as M- 
bands) at 0.75-km spatial resolution and 5 higher-resolution imagery channels (I-bands) at 
0.375-km resolution. With central wavelengths from approximately 0.4 to 12 pm, the 
sensor was designed to be the next-generation global weather and climate imager for 
afternoon polar orbiting observations. To infer cloud properties from the VIIRS 
observations, an accurate and efficient forward radiative transfer model (RTM) is 
invaluable for generating simulated reflectances or brightness temperatures for a variety of 
atmospheric cloud and surface conditions, and can be used for retrieval error analyses and 
instrument calibration efforts. 

Many rigorous radiative transfer schemes, such as the line-by-line radiative transfer 
model (LBLRTM) [Clough et al., 1992, 2005], the adding-doubling (AD) algorithm 
[Twomey et al., 1966; Hovenier, 1969; Hansen, 1971; de Haan et al., 1987], and the 
discrete ordinates radiative transfer (DISORT) method [Chandrasekhar, 1960; Liou, 1973; 
Stamnes et al., 1988; Thomas and Stamnes, 1999], have been developed and widely applied 
under different circumstances. Lor remote sensing problems involving hyperspectrally 
resolved or channel-averaged radiances, a rigorous approach is needed to independently 
and efficiently perform large numbers of simulations for each wavelength/wavenumber 
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due to the significant spectral variations exhibited by molecular absorption, and, in the case 
of channel-averaged simulation, to conduct subsequent spectral integration. Performing 
hundreds or even thousands of monochromatic simulations is extremely time-consuming; 
thus, rigorous RTMs are significantly limited in satellite remote sensing applications 
because of the large number of spatial and temporal observations. Developing 
computationally efficient RTMs for specific satellite-based instruments is critical for 
operational retrievals of atmospheric profiles as well as for the advancement of more 
sophisticated cloud or aerosol property retrievals [Dubuisson et ah, 1996, 2005; Weisz, 
2007; Gamier et ah, 2012, 2013; Chen and Huang, 2014]. 

One approach to achieve computational efficiency is to relax the accuracy constraints 
for each individual simulation; models using two streams [Meador and Weaver, 1980], low 
orders of scatterings [Natraj and Spurr, 2007], or pre-computed look-up-tables [Wang et 
ah, 2012] are typical examples of such approximations. These models normally keep the 
same number of single simulations, but accelerate each individual simulation. Other 
approximate approaches for hyperspectral or band-average calculations such as the 
correlated k-distribution (CKD) method [Arking and Grossman, 1972; Lacis and Oinas, 
1991; Fu and Liou, 1992; Kratz, 1995; Shi et ah, 2009], the principal component method 
[Liu et ah, 2006, 2009], and the optimal spectral sampling method [Moncet et ah, 2008], 
are designed to minimize the computational effort by reducing the number of individual 
radiative transfer simulations within a spectral band. Instead of considering each 
wavenumber within a band, such models conduct rigorous radiative transfer calculations 
at only a small number of representative monochromatic wavenumbers. The two types of 
approaches have been applied to remote sensing, radiative transfer, and GCMs in both the 
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solar and IR spectral regions. In this study, we combine the two approaches to maximize 
computational efficiency. 

This study develops a fast RTM, i.e., an instrument simulator for cloudy atmospheres, 
and uses both of the previously described methodologies to minimize the computational 
time. While the emphasis here is in simulating the VIIRS solar and infrared channels, the 
approach can be easily extended to other imagers. To obtain top-of-atmosphere (TOA) 
reflectances or brightness temperatures, the atmospheric gaseous transmissivity and cloud 
optical properties are essential parameters for a RTM, and will be carefully considered in 
the fast model. Section 2 describes the development of the CKD models to detennine 
gaseous transmissivity, and Section 3 discusses the cloud optical properties. The fast 
radiative transfer models (FRTMs) for VIIRS solar and IR channels, respectively, are 
described in Section 4. Section 5 validates the simulator and compares the simulated results 
with VIIRS observations, and Section 6 summarizes the study. 

2. Determination of gas transmissivity 

As an approximate technique to line-by-line calculations such as those of LBLRTM 
[Clough et ah, 1992, 2005], the CKD is a highly efficient model to account for gaseous 
absorption, and can be easily incorporated into calculations of multiple scattering in 
aerosols and clouds. The CKD model replaces the integral of gas transmissivity over highly 
variable spectral space by one over a much smoother absorption coefficient space. 

The transmissivity of a single gas at constant pressure and temperature within a small 
spectral interval of interest, e.g., the interval of a VIIRS channel, is defined as: 

(i) 
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where k(v) is the gas spectral absorption coefficient at wavenumber v, and u is the gas 
path-length . To accurately obtain the transmissivity, the line-by-line calculation must be 
performed over a very fine wavenumber grid. This calculation requires significant 
computational time because the absorption coefficient is a highly variable function in 
spectral space. However, the transmissivity does not depend on the ordering of the spectral 
absorption coefficient within a given spectral interval, and, thus, instead of integrating over 
the spectral space, Eq. (1) can be expressed as: 

TcM = dv = S^e-^mdk , (2) 

where /(/c) is the normalized probability distribution function for k(v). To obtain /(/c), 
the range of the absorption coefficient k(v) ([k mini k max ]) and the spectral interval are 

k —k 

divided into N and M uniform subintervals, respectively, with width 8k — max mm , and 

8v — We define k 0 — k min , k n — /c n _ 1 + 8k , and v m — (m — 0.5)5u . Thus, the 
probability function /(/r) can be numerically represented by: 

fikd = k i-i < k(v m ) < ki), (3) 

where W is the weighting function equal to unity if its argument condition is satisfied and 
zero otherwise. By defining the cumulative probability function g(k) — f(k')dk\ or 

^min 

in the discrete format, 

g( k i) = -^Zm^SvWikiVn) < k t ) = £}=i f(kj)8k, (4) 

Eq. (1) can be further simplified as: 

T C h(u) = e~ k ^ u dg, (5) 

where g ranges from 0 to 1, and k(g) is a monotonically increasing and smooth function 
of g. Thus, the spectral integration in Eq. (5) can be evaluated with fewer points in g space 
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between zero and unity compared to the number of wavenumber points required for Eq. 
(1). This approximation is known as the ^-distribution method and can be given by: 


T c h(m) = jg 1 e k ^>dg = £f =1 e fc(5i)u A g t . 


( 6 ) 


To extend the ^-distribution method to realistic inhomogeneous atmospheres, the CKD 
method assumes that the ordering of absorption lines with respect to their strengths is the 
same at different temperature and pressure levels, and 


where P(z) and T (z) are the pressure and temperature of the atmospheric layer between zi 
and Z2. 

While the CKD method has been widely used, the spectral k's from multiple gases 
have typically been treated as uncorrelated so that a product rule approximation 
could be used. Edwards and Francis [2000] suggested an approach to treat overlapping 
lines from different gases as a “single gas” by combining the absorption coefficients of 
multiple absorption gases. The equivalent absorption coefficient for a given mixture of 
(N+l) gases is defined as: 


where k t (v) is the absorption coefficient of the z'th gas, and R[ is the density ratio of the z'th 
gas to the reference gas, i.e., the gas with absorption coefficient of k 0 (v). The dominant 
gas can be regarded as the reference gas, and, through Eq. (8), the absorption of (N+l) 
gases is converted into a single-gas case with equivalent absorption coefficient 
K ( v , R 1 , R 2 , ••• R n ). This process can be repeated for a finite number of gas mixture density 


Tchb 0 = ^/ Aw ex P [~ Q k (v.P(z).T(z))pdz] dv 
= fo exp [- f z J k(g,P(z),T(z))pdz ] dg 


( 7 ) 


K(v,R 1 ,R 2 ,-R n ) = k 0 (v) + Zl ikt(v)R u 


( 8 ) 
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ratios for a fixed set of constituent gases, such that K(v) for any given mixture of these 
gases can be found by interpolation. 

To consider an instrument channel’s (i.e., VIIRS) spectral response function (SRF, i.e. 

s(u)), the spectral transmissivity given by Eq. (1) is rewritten as: 

T ch (u) = 1 -S^ v s(v)e- l ‘Wdv, (9) 

where the normalization factor is S — f s(v) dv. To include the SRF in the CKD model, 
the approach given by Edwards and Francis [2000] is used, and Eqs. (3) and (4) are 
modified as: 

fsikd = < k(_v m ) < kd, (10) 

and 

g s (ki) = < k t ) = £}=i f s (kj)8k. (11) 

In Eqs. (10) and (11), the unifonn weighting 8v is changed into s(v m )8v for each spectral 
subinterval. TOA solar spectral irradiance variations in the channel can also be included in 
the normalization factor in a manner similar to the SRFs. 

Fig. 1 illustrates the treatments of overlapping absorption lines, SRF, and solar spectral 

irradiance for the VIIRS M10 channel centered at a wavelength of 1.61 pin. The gas 

spectral absorption coefficients of CO 2 , H 2 O, and CH 4 are plotted in (a), (b), and (c), 
respectively, with (d) showing the effective absorption coefficient of a mixture with mass 
density ratios of H 2 O and CH 4 to CO 2 of 0.1 and 0.001, respectively. Following Eqs. (10) 
and (11), the effective coefficient of the mixture gas is weighted with the corresponding 
SRF (e) and solar spectral irradiance (SI) (f). Fig. 1(g) shows the sorted effective absorption 
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coefficient as a function of g (blue line); results without considering either the SRF or solar 
spectral irradiance is illustrated for comparison (red line). In g space, the absorption 
coefficient becomes a smooth function, and only 4 intervals (shown by the dashed lines in 
(g)) in this channel are used to determine the transmissivity following Eq. (6). Furthermore, 
differences are obvious for results with and without consideration of the SRF, as shown in 
Fig. lg. 

Based on the theories and techniques described, we construct a CKD model for each of 
the VIIRS channels, considering only up to the three most absorptive gases. The parameters 
used to construct the CKD models follow the work done by Ding et al. [2012] for the 
Geostationary Operational Environmental Satellite R (GEOS-R) Advanced Baseline 
Imager (ABI) solar channels. To build the CKD model for each of the VIIRS channels, the 
absorption coefficients of the gases of interest are obtained from the LBLRTM, and the 
molecular absorption line parameters are based on the 2008 edition of the High Resolution 
Transmission (HITRAN) molecular spectroscopic database [Rothman et al., 2008]. The 
absorption coefficients are calculated at 19 pressure levels and 3 temperatures (i.e., 200K, 
260K, and 320K) using the LBLRTM, and each set of absorption coefficients are sorted 
and binned following the same distribution as that of the reference pressure and 
temperature, chosen as 261 hPa and 260 K, respectively. Each of the VIIRS channels is 
divided into 4 to 16 intervals in g-space, the number of which is determined by the degree 
of the gas absorption and the complexity of the overlapping gaseous absorption. 

At a given g and pressure level P 0 , the absorption coefficient at temperature T is given 
by: 

ln[k(g,P 0 ,T)] = a(g,P 0 ) + b(g,P 0 )(T - 260) + c(g,P 0 )(T - 260) 2 , (12) 
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where the coefficients a, b and c are regression coefficients derived from absorption 
coefficients obtained at the three temperature values of 200 K, 260 K, and 320 K for the 
corresponding g values at the 19 pressure levels. In practice, the absorption coefficient at 
an arbitrary temperature and pressure, k(g,P,T), is found first by solving Eq. (12) first, then 
by linear interpolation between values at the two neighboring pressures. 

Table 1 lists parameters of the CKD model for each VIIRS channel, including the 
central wavelength, wavelength range, absorbing gas species considered, and the number 
of g values. Gases considered in this study include H2O, O2, O3, CO2, CH4, and N2O. Up 
to three different absorbers are considered for each channel, and the O3 continuum 

absorption (designated as “cont.” in the table) is included in the IR channels by using the 
equations given by Roberts et al. [1976]. 

The CKD is employed to produce the transmissivity of the atmosphere, and Fig. 2 
compares the band-averaged transmissivity and weighting profiles with results from the 
LBLRTM to illustrate the accuracy of the model for clear sky calculations. Results for four 
VIIRS channels centered at wavelengths 0.555 pm, 1.61 pm, 3.7 pm, and 10.763 pm, are 
shown. In the simulations, the US standard atmospheric profile is divided into 50 layers, 
each being 1 km thick, and the volume mixing ratios of O2, CO2, CH4, and N2O are assumed 
to be uniform with values of 0.21, 3.8xl0' 4 , 1.8xl0' 6 , and 3.2xl0" 7 . From left to right in 
Fig 2, the three columns correspond to the transmissivity in each layer, relative 
transmissivity error, and weighting function profiles of the CKD models. Fig. 2 indicates 
that the relative errors in transmissivity are less than 0.1% for the four VIIRS channels, 
which, from a practical application perspective, are sufficiently accurate for simulating 
TOA reflectance and brightness temperatures under cloudy conditions. 
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3. Determination of cloud optical properties 

For cloudy sky applications, channel-averaged cloud bulk-scattering properties, for 
both liquid and ice phases, are required to determine the absorption, scattering, and 
emission of cloud layers. Here, single-scattering calculations are perfonned at discrete 
sizes and wavelengths, and averaged over the assumed particle size distribution and the 
SRF for each VIIRS channel. For the solar channels, the solar spectral irradiance is also 
considered. For the IR bands, the Planck function representing the thermal IR emission 
from an opaque cloud at 233 K is used. The details for obtaining the channel-averaged 
properties can be found in Baum et al. [2005]. 

For liquid phase clouds, the single-scattering properties are obtained using the Lorenz- 
Mie theory [Mie, 1908]. The cloud droplet size distributions are assumed to be Gamma 
distributions [Hansen and Travis, 1974] with an effective variance of 0.1, and the bulk- 
scattering properties are calculated for effective radii ranging from 2 to 50 pm. 

For ice phase clouds the model choice is critical, because the microphysical and optical 
properties of ice clouds are very sensitive to particle habits, especially the optical properties 
in solar reflectance channels. The ice cloud model used for deriving the MODIS Collection 
6 cloud products is applied in this study. This model consists of a single ice habit, i.e., a 
severely roughened aggregate with eight solid hexagonal columns, that provides excellent 
spectral consistency (i.e., between optical thicknesses inferred from solar reflectance and 
IR channels) and close agreement with the polarization properties from the PARASOL 
observations [Baum et al. 2014]. The column aggregate single-scattering properties are 
obtained from the ice crystal database developed by Yang et al. [2013]. The Gamma size 
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distribution with an effective variance of 0. 1 is assumed for the ice clouds, and the effective 
diameter ranges from 10 to 180 pm in 10 um steps. 

4. VIIRS simulator 

This study considers both the VIIRS solar reflectance and infrared channels. For solar 
channels, cloud multiple scattering plays a much more significant role than the gaseous 
absorption and emission, because absorption is relatively weak at most VIIRS solar 
channels (except the M9 channel with strong water vapor absorption), and emissions at 
these wavelengths are negligible. However, gaseous absorption and emission become as 
important as cloud and aerosol effects (scattering, absorption, and emission) in the infrared 
channels. Thus, the VIIRS solar and infrared channels use different fast radiative transfer 
models (FRTMs). 

The FRTM developed by Wang et al. [2013a] is chosen to calculuate visible through 
shortwave infrared spectral reflectance. This FRTM uses six independent radiative transfer 
equations to approximate the full radiative transfer processes for a combination of cloud, 
aerosol, or molecular layers; the adding-doubling algorithm implemented on a twisted 
icosahedral mesh is incorporated to account for overlapped cloud/aerosol layers. To further 
improve the FRTM efficiency, the bidirectional reflectance and transmittance distribution 
functions (BRDF and BTDF) of cloud layers with different optical thic kn esses and 
effective particle sizes are pre-calculated using the 128-stream DISORT, and the optical 
properties of water and ice clouds, discussed in the previous section, are included. 
Generally, the FRTM is approximately two orders of magnitude faster than the 128-stream 
DISORT, and obtains TOA reflectance with relative errors normally less than 5%. 
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For the infrared channels, the simulator uses the FRTM developed by Wang et al. 
[2011, 2013b] to obtain the TOA brightness temperature. Similar to the FRTM used for the 
solar channels, the effects of cloud layers are efficiently considered by pre-calculated look- 
up-tables at various optical thicknesses and effective particle sizes, which include the 
reflectance, transmittance, emissivity, and effective temperature. The CKD model 
discussed in Section 2 is used to account for atmospheric gas absorption. The TOA 
brightness temperature differences (BTDs) given by the fast model and the rigorous 
DISORT are less than 0.15 K. Furthermore, the FRTM is more than three orders of 
magnitude faster than the corresponding a DISORT implementation with 32 streams. 

The most current DISORT code (DISORT 2.0 beta) is used to calculate the look-up- 
tables for both FRTMS based on the band-averaged optical properties. In DISORT, the 
cloud phase functions are defined in terms of their Legendre polynomial expansion 
coefficients and, for large ice particles, thousands of Legendre polynomial terms are 
needed to represent a phase function that has a strong forward peak. The number of 
expansion terms can be greatly reduced, however, by truncating the forward peak of the 
phase function. The 5-fit method [Hu et al., 2000] is used here for the solar channels, such 
that only 128 terms are sufficient to reconstruct the phase function, significantly decreasing 
computational time in the DISORT simulations. Corresponding scaling adjustments are 
made to the optical thickness and single-scattering albedo to account for the truncated 
forward energy. 

The solar and IR FRTMs are combined with the CKD method to maximize the 
efficiency of the simulator. For each simulation, the FRTM is performed with the 
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transmissivity obtained under each g value, and the TO A radiance of a given channel is 
obtained by: 

I ch, simulator ~ Xi=i i[T(giMgu (13) 

where T(g t ) refers to the transmissivity of the z'th g value, and I\T(gi)\ is the 
corresponding TOA radiance given by the FRTM. 

The accuracy and efficiency of the simulator, i.e., the model based on the combination 
of the CKD and FRTMs, are evaluated by comparing the simulated band-averaged 
reflectance or BT at the TOA with the rigorous solutions given by the combination of the 
LBLRTM and DISORT. The spectral resolution of the LBLRTM+DISORT simulation is 
chosen to be 1 .0 and 0. 1 cm' 1 for the solar and infrared channels, respectively, and the TOA 
upwelling radiances are then averaged over the spectrum considering the SRF: 


j _ [av s(v)I(v)dv 

*ch,std ~ c > 


(14) 


where / (v) is the radiation at wavenumber v given by the DISORT. With the TOA 
radiance, obtaining the corresponding reflectivity and brightness temperature is 
straightforward. 

A set of comparisons at solar channels is implemented between the fast simulator and 
the rigorous approach, with the relative errors of TOA reflectance as a function of the 
viewing zenith angle shown in Fig. 3. Three VIIRS channels (M4, M10, and Mil) are 

considered for the comparison, and three solar zenith angles values (10°, 30° and 50°) are 


used. The left panels are for an optical thickness of 5, and the right panels are for a value 
of 20. For all cases, the simulator yields relative errors less than 1.5% with respect to the 
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rigorous approach. The relative errors slightly increase as the solar zenith angle increases, 
and show little dependence on the cloud optical thickness. 

To validate the simulator at the IR channels, Fig. 4 illustrates the brightness temperature 
differences (BTDs) given by the simulator and the LBLRTM+DISORT at three IR 
channels (Ml 4, Ml 5, and Ml 6). The BTD is defined as: 

BTD = BTsimulator- BTLBLRTM+DISORT. (15) 

Each panel of Fig. 4 shows the comparison of simulation results based on a surface 

emissivity of 1 and viewing zenith angle of 20° under different cloud conditions. The BTD 

is expressed as a function of optical thickness, and cloud particle effective diameters and 
top temperatures used for the simulations are listed in the figure. The errors in the BTDs 
are smaller than 0.2 K and decrease to less than 0. 1 K for optically thick clouds. 

5. Comparison with VIIRS observations 

We developed a fast radiance simulator to calculate TOA reflectances or brightness 
temperatures of a cloudy atmosphere based on a combination of the CKD and FRTMs for 
the VIIRS solar and infrared channels. This section highlights a case study to assess the 
performance of the simulator by comparing simulated TOA reflectances and brightness 
temperatures with those from VIIRS observations. 

As discussed in the previous sections, the forward model requires atmospheric profiles 
and cloud properties as input parameters. For this comparison, atmospheric profile data is 
obtained from the Modem Era Retrospective-analysis for Research and Applications 
(MERRA) [Rienecker et ah, 2008] instantaneous 3-hour vertical atmospheric profile 
product (i.e., inst3_3d_asm_Cp) that provides temperature, geopotential height, water 
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vapor, and ozone concentrations at 42 pressure levels on a 288 x 144 mesh grid with a 1 .25 


x 1 .25° resolution. The cloud input to the simulator, including cloud thermodynamic phase, 

top pressure, optical thickness, and effective particle size, is from the operational MODIS 
Collection 6 cloud products (i.e., the Aqua MODIS MYD06 L2 product), and the 1 km- 
resolution geolocation is obtained from the MYD03 dataset. The atmospheric profiles and 
cloud properties are collocated with the VIIRS observations, with VIIRS solar and sensor 
view geometries used as simulator inputs. 

A flowchart summarizing the fast VIIRS radiance simulator process is shown in Fig. 
5. The atmospheric profile is input into the CKD models to generate transmissivity of 
absorption gases, and also provides the temperature profile for the IR simulator. The 
gaseous transmissivity and cloud properties are used by the FRTMs to calculate 
reflectances or brightness temperatures of given VIIRS solar-satellite geometries, and the 
simulated results are compared with the corresponding VIIRS measurements to assess the 
performance of the simulator. 

To avoid uncertainties associated with satellite-based cloud retrievals due to larger 
surface reflectances over land, only observations over ocean are considered for the present 
case study, specifically, the VIIRS granules taken over the South Pacific Ocean from 
00:38:47.4 to 00:41:38.2 UTC on 3 January 2014. The collocated MODIS granule is taken 
at 00:45 UTC, approximately 4 to 6 minutes behind the VIIRS observations. Fig. 6 shows 
the true color images of the (a) MODIS and (b) VIIRS granules, with the MODIS cloud 
optical thickness and effective particle radius shown in (c) and (d), respectively. The red 
boxes in the RGBs show the region of the simulation, which is largely covered by ice 
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clouds. In the simulation region, the cloud optical thickness ranges from a few to over 50, 
and the effective particle radius shows values from 5 to approximately 40 pm. 

Fig. 7 compares the observed (left panels) and simulated (right panels) reflectances at 
VIIRS 0.86- and 2.25-pm channels (M7 and Mil). The reflectance values at the 0.86-pm 
channel are larger than those at the 2.25-pm channel, which is mainly due to the significant 
differences in the scattering properties of ice clouds at the two channels (e.g., ice clouds 
are much more absorptive at the 2.25-pm channel). The patterns of the reflectances given 
by the fast simulator are almost the same as the VIIRS observations at each channel. 
Furthermore, the agreement indicates the performance of the MOD06 product that is used 
as the input parameters. However, note that the liquid water and ice cloud radiative models 
used in the FRTM were specifically chosen to match the models used in MYD06 L2 (see 
Sect. 3). 

Fig. 8 is the same as Fig. 7 but for brightness temperatures at three VIIRS IR channels 
(i.e., 8.55-, 10. 76-, and 12.01-pm channels (M14, M15, and M16) from upper to lower 
panels). The simulated brightness temperatures show close agreement with the 
observations. However, noticeable differences exist in some regions of the granules, and 
this may be due to the uncertainties in either atmospheric profiles or cloud height. The ice 
cloud properties, which are retrieved from MODIS solar-channel observations, may also 
yield significant errors when applied to IR channels, because the cloud optical thic kn ess 
and effective particle radius inferred from the solar and IR channels can be quite different 
[Baum et ah, 2014]. The case study indicates the exceptional performance of the VIIRS 
simulator for both solar and IR channels. 
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6. Summary 

This study developed a computaionally efficient simulator for the VIIRS solar and IR 
channels in cloudy atmospheres that can be used in cloud property evaluations and 
retrievals. The absorption of atmospheric gases and overlapping gas absorption is 
accounted for using a VIlRS-specific CKD that considers both the spectral response 
function and solar spectral irradiance. The accuracy of the transmissivity profile is 
estimated by comparing with the exact line-by-line results, and the relative errors in 
transmissivity are less than 0.1% for all VIIRS channels. Two fast RTMs are used to 
consider absorption, scattering, and emission of cloud layers. The channel-averaged bulk- 
scattering properties of roughened hexagonal columns are used for ice cloud, and the 
properties of water cloud are given by the Lorenz-Mie theory. By comparing with the 
rigorous DISORT results, the relative errors for TOA reflectance at VIIRS solar channels 
are less than 1.5%, and the differences in brightness temperatures at the IR channels are 
less than 0.25K. The present simulator is more computationally efficient than the standard 
LBLRTM+DISORT by over three orders of magnitude. With collocated MERRA 
atmospheric profiles and cloud optical thickness and effective particle diameter from the 
MODIS cloud product as input, the reflectances and brightness temperatures calculated by 
the fast simulator show close agreement with concurrent VIIRS solar and IR observations. 
Considering the accuracy and efficiency provided, the simulator can be used directly for 
cloud property retrievals related to VIIRS observations. While our fast radiative transfer 
model (FRTM) used in this study was applied to VIIRS channels, the FRTM can also be 
developed for MODIS and other satellite or airborne imagers with similar spectral 
coverage. 
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534 TABLE 1. CKD model parameters for the VIIRS channels. 


VIIRS 

Channels 

Central 

Wavelength 

(pm) 

Wavelength 
Range (pm) 

Gas Considered 

Number of 
g values 

Ml 

0.412 

0.402-0.422 

H 2 0 

4 

M2 

0.445 

0.436-0.454 

H 2 O 

4 

M3 

0.488 

0.478-0.488 

H 2 O 

16 

M4 

0.555 

0.545-0.565 

H 2 O 

4 

M5 

0.672 

0.662-0.682 

H 2 O, O 2 , O 3 

8 

M6 

0.746 

0.739-0.754 

H 2 O, O 2 

8 

M7 

0.865 

0.846-0.885 

H 2 O, O 2 

8 

M8 

1.240 

1.23-1.24 

H 2 O, O 2 , CO 2 

8 

M9 

1.378 

1.371-1.386 

H 2 O 

12 

M10 

1.61 

1.58-1.64 

H 2 O, CO 2 , CH4 

4 

Mil 

2.25 

2.23-2.28 

CH 4 , H 2 O, N 2 O 

4 

M12 

3.7 

3.61-3.79 

H 2 O, CH 4 , 03 

4 

M13 

4.05 

3.97-4.13 

C02, H 20 , N 20 

16 

M14 

8.55 

8.4-8. 7 

H 2 O, N 2 O, O 3 , cont. 

8 

M15 

10.763 

10.26-11.26 

H 2 O, CO 2 , O 3 , cont. 

4 

M16 

12.013 

11.54-12.49 

H 2 O, CO 2 , O 3 , cont. 

8 

11 

0.64 

0.60-0.68 

H 2 O, O 2 

4 

12 

0.865 

0.85-0.88 

H 2 O, O 2 

8 

13 

1.61 

1.58-1.64 

H 2 O, CO 2 , CH 4 

4 

14 

3.74 

3.55-3.93 

H 2 O, N 2 O, CH 4 

4 

15 

11.45 

10.5-12.4 

H 2 O, CO 2 , O 3 , cont. 

16 
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539 FIG. 1 . Example of the treatment of overlap absorption lines, spectral response function 


540 (SRF), and solar irradiance (SI) based on the CKD algorithm for the VIIRS 1.61 -pm 


541 channel, (a)-(c) Absorption coefficient as a function of wavenumber for H 2 O, CO 2 , and 


542 CH4..(d) Absorption coefficient for a mixed case with the ratio of the three gases being 


543 1:0.1:0.001. (e) Spectral response function, (f) Solar irradiance at the top of atmosphere. 


544 (g) Absorption coefficient as a function of cumulative probability for the mixed gas with 


545 and without the inclusion of the spectral response function and solar irradiance. 
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FIG. 2. Band-averaged transmissivity (left) calculated from the LBLRTM and CKD 
models for VIIRS 0.555-, 1.61-, 3.7-, and 10.763-pm channels. Corresponding relative 
errors (middle) and weight function (right) profiles of the CKD models. The U.S. standard 
atmospheric profile is used in the calculations. 
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FIG. 3. Relative errors of TOA reflectance at M4 (0.555 jam), M10 (1.61 jam), and M12 
(3.7 jam) channels by the VIIRS simulator in comparison with the benchmark model. 
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FIG. 4. Brightness temperature difference (CKD+FRTM - LBL+DISORT) as a function 
of ice optical thickness for a viewing zenith angle of 20° at three VIIRS infrared channels 
(8.55-, 10. 763-, and 12.013-pm). 
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561 FIG. 5. Flowchart outlining the fast VIIRS radiance simulator using the collocated 


562 MERRA atmospheric profile and MODIS retrieved ice cloud thic kn ess and effective 


563 particle size. 
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(a) MOD1S RGB Image 
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(b) VIIRS RGB Image 
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FIG. 6. (a) MODIS RGB image for a granule on 3 January 2014 at 0:45 UTC taken over 


(c) MODIS C6 COT 
-170 -160 -150 


(d) MODIS C6 Effective Radius (gm) 
-170 -160 -150 


564 

565 


566 the South Pacific Ocean, (b) RGB image of the corresponding VIIRS granules from 


567 00:38:47.4 to 00:41:38.2 UTC of the same day. (c) and (d) are the MODIS Collection 6 


568 cloud optical thickness and particle effective radius of ice clouds. 
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(a) Observed M7 Reflectance 
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(b) Simulated M7 Reflectance 
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(d) Simulated Mil Reflectance 
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FIG. 7. Comparison between the observed (left panels) and simulated (right panels) 
reflectances at the VIIRS M7 (0.865 pm) and Mil (2.25 pm) channels. 
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(b). Simulated M14 BT (K) 


-170 -165 -160 -155 

(d). Simulated Ml 5 BT (K) 


-170 -165 -160 -155 
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573 

574 FIG. 8. Comparison between the observed (left panels) and simulated (right panels) 

575 brightness temperatures at the VIIRS M14 (8.55 pm), M15 (10.763 pm), and Mil (12.013 

576 pm) channels. 
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